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CN| ' Abstract. Together with the concept of reversibility, another relevant physical 

, notion is time-symmetry, which expresses that there is no way of distinguishing 

between backward and forward time directions. This notion, found in physical 
£^ | theories, has been neglected in the area of discrete dynamical systems. Here we 

formalize it in the context of cellular automata and establish some basic facts and 
relations. We also state some open problems that may encourage further research 
on the topic. 
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1. Introduction 

An important property that may be present or not in physical or abstract dy- 
namical systems is reversibility; consequently, it has also been an active topic of 
research in the context of cellular automata[S]. At least two particular reasons for 
' this interest are often mentioned: on one hand, if CA are seen as models for massive 

cn ■ distributed computation, then Landauer's principle suggests that we should focus on 

reversible cases. On the other hand, reversibility is often observed in real systems; 
it is therefore desirable in models of them[T5]. Furthermore, a number of interesting 
results (like the dimension- sensitive difficulty of deciding reversibility [6]) have kept 
reversible CA in sight over the years. 

However, there is one aspect of reversibility, as seen in real systems, which 
has been mostly neglected when considering cellular automata (in fact, for discrete 
dynamics in general): the dynamical laws governing physical reality seem to be 
not only reversible, but time- symmetric. For Newtonian mechanics, relativity or 
quantum mechanics, we can go back in time by applying the same dynamics, pro- 
vided that we change the sense of time's arrow, through a specific transformation 
of phase-space. In the simplest example, Newtonian mechanics, the transformation 
leaves masses and positions unchanged, but reverses the sign of momenta. 
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In the most general sense, we say that a dynamical system (X, T) is time- 
symmetric if there exists a reversible R : X — > X such that R o T o R^ 1 = 
T _1 [10J (notice that this applies to systems with discrete or continuous time). How- 
ever, time-symmetries observed in physical systems follow usually a more restricted 
definition, in which R~ l = R, and therefore R is an involution on X. This is a 
natural restriction, which follows whenever there is no way to distinguish where the 
arrow of time is heading. Apparent irreversibility (Loschmidt's paradox) comes only 
from macroscopic (i.e., coarse-grained) differences in entropy. 

Here we will discuss some basic facts about time-symmetric cellular automata, 
defined as those CA F for which there exists an involution H (which is a CA itself) 
such that 

F- 1 = HoFoH (1.1) 

Requiring if to be a CA is somewhat arbitrary, since for other systems the time- 
reversing transformation is not necessarily of the same nature as the dynamics (in 
fact, the physical theories discussed above are continuous in time). The reason for 
this restriction is that we expect reversibility (including the particular case of time- 
symmetry) to be a local property. Even if we do not address the case when H is not 
a CA, it may be an interesting direction for future studies. 

CA are usually defined over a full shift S z , but they can also be studied over 
(stable) subshifts. We remark that in this case, for time-symmetry to apply, the 
subshift must be stable for both F and H. This may cause some problems, since 
subsystems of a time-symmetric CA cannot be assumed to be time-symmetric too, 
even if they are stable for F. 

2. Some motivating examples 

Not only our models of physical reality turn out to exhibit time-symmetry; it 
is also found in some well known reversible discrete dynamical system. We show in 
this section how it applies to two 2D system, Margolus' billiard and Langton's ant. 
As a technical note, notice that in both cases the system is not originally described 
as a cellular automaton in the strict sense; therefore, we describe for each of them 
a CA that contains them as particular case for a subshift of valid configurations. 
This should -in principle- be followed by an extension of the rule to the full shift, 
and such that the system remains time-symmetric; however, doing that is not really 
required, since we want to show the time-symmetry of the original system; we hence 
restrict ourselves to the valid subshift. 

Margolus' billiard. A well known example of time-symmetric CA is the 
Billiard ball model of Margolus [12]. It is not a proper CA, but rather a so-called 
partitioned CA, where the space Z 2 is partitioned in 2 x 2 blocks of cells in two 
different ways (see Figure QJa)). A transformation is applied to each block of each 
partition alternately. It is easy to see that such an automaton is reversible if and 
only if its local transformation is one-to-one. The rule used by Margolus is shown in 
Figure [T|b). It tries to emulate balls that move in straight lines, colliding elastically 
with each other or with static obstacles. The importance of this model comes from 
its Turing- universality, proved in [12] by computing reversible Fredkin gates [3J. 

We can express Margolus' system in terms of a CA with alphabet {white, black] x 
{/\\, v /,\} and Moore neighbourhood. Here the first layer (the white/black 
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component) represents the states of the original Margolus model, and the other rep- 
resents the current partition, along with the relative position of the cell within its 
current block. This layer must be initialized in an appropriate way in order to work 
correctly (see Figure H])(c). 

Notice that reversing the arrows makes the partition flip to the alternative one. 
At each time step, each cell computes its next white/black state by applying Mar- 
golus' rule to the quadrant indicated by its arrow, and then reverses its arrow. Each 
of this actions -the first on the first layer, the second on the second- is an involution. 
Furthermore, if at one time step we omit any of them, further iterations will make 
the automaton evolve back in time. 
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Figure 1: (a) The two partitions of the Margolus model are shown, one with solid 
lines and the other with dashed lines, (b) The Billiard Ball Model is 
defined through a permutation over 2x2 blocks of cells, (c) The current 
partition is obtained by grouping the four cells that point to the same 
point; reversing the arrows gives the alternate partition. 



Langton's ant. Langton's ant was introduced in [TTJ together with several 
models emulating different life properties. It was also defined in physics as a model 
for particles presenting self correlated trajectories pQ. The model can be seen as a 
Turing machine working on a 2-dimensional tape. Its internal state is an arrow that 
represents its last movement direction. At each step, the ant turns to the left or to 
the right depending on the cell color (white or black), it flips this color and moves one 
cell forward (see Figure [21(a) ) . Besides being Turing- universal J2], its celebrity is due 
mostly to its particular behavior over finite initial configurations. Simulations show 
that it always falls eventually into a repetitive movement -of period 104- that makes 
it propagate unboundedly (see Figure 12(b)); this assertion has not been proved, and 
appears to be very difficult despite the simplicity of the transition rule. 

Langton's ant can also be described in terms of a CA with Moore neighborhood 
and state cell {head, tail, empty} x {white, black}. We represent the arrow through 
two adjacent cells, one in state head and the other in state tail. The cell in state 
tail always becomes empty, while the cell in state head always becomes tail and flip 
its color. Cells adjacent to a head can decide to become head themselves by looking 
at the tail position and the color of the head cell. The system simulates Langton's 
ant only if it starts with only one ant. 
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Here again, we can define the involution consisting in exchanging tails and heads. 
This immediately makes the ant come back to the cell it just had left, which it finds 
in the color opposite to the one it had found before, causing the ant to turn in the 
opposite direction, which in turn makes it again go to a previously visited cell, and 
so on: the ant will forever retrace (and undo) its past trajectory. 



5] 




Figure 2: (a) Langton's ant rule, (b) Space configuration at iteration 10,837 after 
starting with every cell in white color. 



3. Basic results 

Proposition 3.1. Let F be a CA. Then the following are equivalent: 

(1) F is time-symmetric. 

(2) There exists an involution H such that (F o H) is an involution. 

(3) F is the composition of two involutions. 

Proof. (1) (2) 



Let F and H be the CA satisfying 11.11 Then 

(F0H) 2 = F0H0F0H = FoT 1 = id 

(2) => (3) 

Take H from (2) and let G = F o H which is an involution. We have 
F = Foid = Fo(HoH) = (FoH)oH = GoH 

(3) =► (1) 

Let G and H be involutions such that F — G o H . Then 

F- 1 = (GoH)- 1 = H^oG- 1 = HoG = H0G0H0H = H0F0H 



184 



A. MOREIRA AND A. GAJARDO 



Remarks 3.2. The following additional facts are noteworthy: 

(1) If F is time-symmetric, then so is its inverse F _1 . Moreover, if F = Go H is 
a decomposition into involutions, then F^ 1 = H o G is a decomposition for 
the inverse. If H was the involution verifying 11.11 then G plays that role for 
F- 1 . 

(2) For any i G Z, F l is also time-symmetric. 

(3) The identity is a (trivial) involution; from there and the third condition we 
have that any involution is trivially time-symmetric. 

(4) Not every reversible CA is time-symmetric. For example, a (the shift): if for 
some H, (a o H) o (<j o H) = id, since any CA commutes with the shift, we 
would have a 1 = id, which is a contradiction. 

The following diagram commutes: 



X 



X 



F 



-i 



X 



X 



Moreover, if we use F = G o H to decompose the dynamics into the alternate 
applications of the involutions, so that successive configurations are computed as 

, where both F and 
This curious situation is 



c' t = H(c t ), c t +i = G(c' t ), we get a dynamics cq, 



F 1 are being iterated: Ct+i 
represented in Figure El 



'0) 



F(ct) and c' t 



Ci, c 



i- 



t+i 




Figure 3: The decomposition of time-symmetric CA into alternating involutions cre- 
ates a situation where both F and its inverse can be read from the space- 
time diagram as time moves forward (or backward). 
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4. Involutions 

Involutions are quite infrequent in the space of CA. For example, the only ele- 
mentary CA of period two are the identity and its negation. 

It is easy to decide whether a given CA is an involution or not: we just have 
to compute its square and compare it to the identity. Nevertheless, if we want 
to enumerate the set of involutions, this procedure is very slow. A constructive 
characterization or a practical set of strong necessary conditions is still missing. 

Meanwhile, it may be useful to consider restricted families of CA. For instance, 
if we restrict ourselves to additive CA, we get an alternative characterization in 
terms of coefficients. If we consider an additive CA of radius r defined on (Z m ) z by 
the local rule 



h{x— r , . . . , Xq, . . . , x r ) ^ 



CljXj 



then by applying the rule twice and grouping the terms it is easily seen that h is an 
involution if and only if 

r r+min{j,0} 

a + 2 flia_j = 1 and cuaj-i = , Wj ^ 

i=l i=— r+max{j,0} 

For instance, for m = 4 and r = 1, we get 2a_i<2o = 2a ai = = a\ = and 
a,Q + 2a_iai = 1. Putting a_i = a± = 2 we obtain all the zeroes, and with ao = 1 or 
ao = 3 we have the last condition. An example of an additive involution is thus 

h(x^i, x , xi) = (2(x_x + x\) + 3x ) mod 4 (4.1) 

Another well studied family of CAs are permutative ones. Unfortunately, we do 
not have a characterization there. A necessary condition is given by the following 
fact: 

Proposition 4.1. Given an involution h, the following two assertions are equiva- 
lent: 

• h is left-permutative. 

• h is one-way to the righ^\. 

Proof. If h is left-permutative of left radius I, h 2 is also left-permutative of left 
radius 21, but the left radius of h 2 is 0, then / is 0. Conversely, if h is oneway to the 
right and it is not permutative, there exists x\...x n and y\ such that h(xi, ■■■,x n ) = 
h(yi,x 2 , --,x n ). Taking any extension z of x\..x n to Z, and z[ = Zi for every i ^ 1 
and z[ = yi, we have that h(z)[i ;00 [ = h(z')[i t00 \. Thus, z[ = h 2 (z) 1 = h 2 (z')i = z±, 
which is a contradiction. ■ 

Thus the involution in ( 14. ip . which is clearly two-way, is not permutative. An 
example of a permutative involution of radius r is: 

{1 if x_ r = A x = 2 
2 if x- r = A x = 1 (4.2) 
xq otherwise 

This kind of construction is the simplest one, with permutations which are trans- 
positions and which do not affect the states that regulate their application. More 



1 The neighbourhood is a finite subset N of Nq. 



186 



A. MOREIRA AND A. GAJARDO 



complicated examples may have associated permutations which are not transposi- 
tions. 

Example ( I4.2p . incidentally, shows how an involution can have an arbitrarily 
large neighbourhood. This long-distance dependence may be lost when it is com- 
posed with another involution, making the determination of time-symmetry non- 
trivial. For instance, the involution (I4.2p yields the permutation (12) of radius 
when it is composed with 

f 1 if x_ r ^ A x = 2 
g(x_ r , ..x ) = < 2 if x_ r 7^ A x = 1 
[ xq otherwise 

5. Diversity in the class of time-symmetric CA 

One simple example of time-symmetric automata is given by the following. 

Proposition 5.1. Every reversible CA of radius is time- symmetric. 

Proof. Let / be the local rule of a reversible CA of radius 0, and let S be its set 
of states. Suppose first that / : S — > S is a cyclic permutation and, without loss 
of generality, that S = {0, ..,n — 1} and f(i) — % + 1 mod n. Let us define the 
involution h(i) = n — i — 1. Consider g = h o f; for % < n — 1, g(i) = h(f(i)) = 
h(i + 1) = n — (i + 1) — 1 = n — i — 2, while otherwise g{n — 1) = h(f(n — 1)) = n — 1. 
Thus g is an involution: g 2 (n — 1) = g(n — 1) = n — 1, and for other i, g 2 (i) = 
g(n — i — 2) — n — (n — % — 2) — 2 = i. Since / = h o h o f=h o g is the composition 
of two involutions, it is time-symmetric. 

If / decomposes into more than one cycle, we define h and g as before over each 
of them, obtaining again a decomposition into involutions. ■ 

It is important to notice here the preservation of time-symmetry under conju- 
gacy. 

Proposition 5.2. If F is conjugated to T andT is time-symmetric, then F is also 
time-symmetric. 

Proof. From time-symmetry, there is an involution H such that T _1 = H o T o H. 
From conjugacy, there is a bijective, continuous, shift- commuting cj) such that T = 
cj) o F o <p~ l . Then we have (removing the composition symbol, for clarity) that 

F- 1 = cfy^T- 1 ^ = (j)~ x HTH(f) = ^H^F^Hcj) = GFG 

and G = <f>~ x H§ is cleary an involution, making F time-symmetric. ■ 

Periodic CA of radius r > behave almost like a CA with radius 0, in the 
sense that information cannot travel "very far"; this makes them nearly time- 
symmetric, because they are conjugated to a subshift of a radius CA. To see 
this, let F be a p-periodic CA with states S and define <p : S z — > (S p ) z as 
ip(x)i = (xi, F(x)i, .., F p ~ 1 (x)i). This tp is continuous and injective, and the in- 
duced CA F' in (S p ) z has radius and period p; its local rule is /'(ao, ai, •-, a p -i) = 
(ai, a 2 , .., a p _i, a ). From Proposition 15.11 we see that F' is time-symmetric; more- 
over, <p(S z ) is .F' invariant. However, f(S z ) is not invariant for the involution 
defined in the proof, and therefore we cannot conclude that F is time-symmetric. 
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When we consider the group of reversible CA with the composition operation, 
involutions correspond to elements of order two. These objects were already con- 
sidered by Hedlund et aim [5], where the last theorem shows that the composition 
of two involutions (i.e., time-symmetric CA) can have infinite order. The example 
that proves this theorem is defined on alphabet {0,1} and consists in the compo- 
sition of a, the negation of the identity, and /3, a CA that negates Xi if and only 
if Xj_iXj + iXj + 2 = 101. Figure H] shows simulations of this CA; the Hedlund's proof 
exhibited a configuration with infinite orbit consisting in one traveling signal over a 
periodic background like the one appearing in Figure 0Ja). 
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Figure 4: Simulations of the time-symmetric CA defined by Hedlund et al on pe- 
riodic boundary conditions, (a) Two traveling signals, (b) A simulation 
over a random initial configuration. 



This already suggests a variety of dynamical behaviors within the class. But 
the examples given in Section |2] are even more interesting, as they correspond to 
Turing-complete systems. The following results shows that, indeed, the whole range 
of reversible dynamical behaviors can be observed in time-symmetric CA. 

Proposition 5.3. Let F be a ID reversible CA. Then there exists a ID CA F which 
is time- symmetric and simulates F in real time. 

Proof. Let / be the local rule of F and denote with f~ l the local rule of its inverse 
F^ 1 ; let £ and r be large enough so that iV = {—£, . . . ,r} contains the neighbour- 
hoods of both / and f~ l ; finally, let S be the set of states. We define the CA F 
with neighbourhood iV and states S 2 , through the local rule 

f((x- t ,y-t),...,(x r ,y r )) = (f(x-e, ■■■ ,x r ),/ _1 (^, ... ,y r )) 

F simulates F in real time: to project the space-time diagram of F into that of F, 
we just discard the second component of the ordered pairs. By discarding the first 
component instead, we note that F simulates F" 1 as well. 

Let H be the involution given by the radius local rule h(x, y) = (y, x). Abusing 
notation, denote configurations c G (S 2 ) z as pairs (x,y) 6 (S z ) 2 . Then we have 

FoH(x,y) = F(y,x) = (F(y), F~\x)) 

and 

(FoH) 2 (x,y) = F o H(F(y),F~ 1 (x)) = F(F-\x), F(y)) = (x,y) 
and thus F is time-symmetric. ■ 
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Cellular automata are said to be intrinsically universal if they are able to simu- 
late any other CA. The details vary according to the accepted notion of simulation, 
from which there is a variety. Delorme et al [2] have recently reviewed and com- 
pleted the study of three of these, surjective, injective and mixed simulation, and 
shown that for every pair of CA F and G, F x G simulates both F and G in all 
three senses. 

Corollary 5.4. There exist time-symmetric CA which are intrinsically universal 
within the class of reversible CA. 

Proof. This follows from the previous results and comment, and from the existence 
of reversible intrinsically universal CA (see for example [33]). ■ 

Notice that reversible CA cannot simulate arbitrary CA: intrinsic universality 
is therefore limited to the reversible class, and time-symmetric CA are as general as 
reversible CA can get. Turing- universality is not limited by reversibility (information 
can be "swept away" to preserve it and maintain reversibility) and hence is implied 
by reversible intrinsic universality. 

Not every reversible CA is time-symmetric; a simple example is the shift a, 
which commutes with every CA and therefore cannot satisfy equation fll.ip for any 
involution H. But in general, it is not easy to prove non-time-symmetry. An in- 
teresting theory which may provide better tools for doing this, and possibly for 
characterizing time-symmetry, is the one developed by Kari in [?]. We will not re- 
produce here his construction, but one important fact is the following: he introduces 
a morphism h_ from the set of reversible CA with the composition (Aut(A), o) into 
the set of rational numbers with the multiplication (Q, •): h_(f o g) — h_(f)h_(g) 
for all reversible /, g. Clearly, involutions and every periodic CA are in the kernel of 
h_. Moreover, since time-symmetric CA are compositions of involutions, they are 
in this kernel as well. We do not presently know whether they are identical to the 
whole kernel or not. 

Kari proves that reversible CA which are not in this kernel are compositions of 
some element of the kernel with a partial shift which is easily computed from the 
value of h_; in turn, every element of the kernel can be written as a composition of 
two block permutations (akin to Margolus rule), and thus he expresses reversible CA 
in an explicitly reversible way. Our motivation here is different, but the approach is 
promising and the connection should be explored. As a first conclusion, we obtain 
that every CA in the kernel of h_ is a composition of two time-symmetric CAs, and 
hence is also the composition of four involutions. 

6. Conclusions 

We believe that this note just scratches the surface of the topic of time-symmetry 
in cellular automata; their rich internal structure and the connection to physical 
models suggests that much more can be done with them. 

On the other hand, as shown by the examples in Section [21 time-symmetric 
CA are actually quite familiar to CA researchers, and have appeared in different 
contexts. Some cases are very explicit, like the automata constructed in the proof of 
undecidability of periodicity [5], which actually include an "arrow of time" toggle. 
Moreover, there are ways of constructing CA rules that make the construction of 
time-symmetric CA straightforward. For instance, Margolus' billiard is an example 
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of a block automata, i.e., a system which is a composition of two functions applied 
to independent blocks of the configuration. By incorporating the current function 
and block to be applied into the configuration, a block automaton can always be 
expressed as a CA. Defining an involution that toggles the current block is a good 
idea to prove time-symmetry, but it only works if both block functions are also 
involutions. What must be stressed is that this is only a sufficient condition; the 
system may be time-symmetric by means of an entirely different involution. 

Likewise, partitioned CA (in the sense of Morita [13]) can easily give birth to 
time-symmetric CA. In that case, cells are partitioned into sub-cells, one for each 
neighbours; iteration proceeds by the alternation between an exchange step, where 
cells exchange the contents of the sub-cells associated to each other, and a step 
which applies a block transformation on the cell. This scheme was succesfully used 
to construct reversible CA (all we need is a reversible block transformation), and 
can produce time-symmetric CA as well if the block transformation is chosen as 
an involution: the exchange step already is one. Again, what we want to stress is 
that this is a sufficient condition: we could have a partitioned CA which is time- 
symmetric while having a non-involutive block transformation, if the decomposition 
happens to be another one. 

There are several interesting questions that should probably be addressed next, 
and have appeared along this text: 

• Is there a constructive characterization of CA involutions that can make 
their enumeration practical? Right now the only way we have to find the 
involutions is to test all CA exhaustively; some trivial necessary conditions 
can be used to reduce the search, but they are not enough to make it efficient. 

• Is time-symmetry a decidable property? Since the definition calls for the 
existence of an involution that verifies a condition, a bound on the necessary 
neighbourhood for the involution would be enough to ensure decidability. 

• Do time-symmetric CA correspond to the kernel of Kari's h- morphism? 
We conjecture a positive answer for these three questions, at least in dimension 

1; Kari's result on undecidability of reversibility in dimension 2 [6] suggests that 
answers here may be dimension- sensitive too. The answers to the questions may 
be related to each other. For instance, a better understanding of the structure 
of involutions may be useful for bounding the required neighbourhood and thus 
deciding time-symmetry. On the other hand, since h_ is easily computed, a positive 
answer to the third question would imply a positive answer to the second as well. 

Notice that if the answer to the third question is positive, then time-symmetric 
CA would be closed under composition. This is by no means obvious, and in fact it 
is a further interesting open question. 

A further direction for future work may be the study of time-symmetry in other 
discrete dynamical systems. In each case an important issue is to precise what kind 
of involution is to be applied. Generally speaking, what we need is an involutive and 
hopefully local transformation of the system's configuration. That transformation 
may not be, in general, an object of the same kind as the dynamics itself: that was 
the case for CA because of the special nature of CA, which transform the whole 
configuration in discrete time too, and will be the case for automata networks in 
general. In other cases, like for instance Turing machines, it is not only difficult (the 
composition of two Turing machines moves the head two steps, and is no longer a 
Turing machine unless we extend the definitions) but also not expected; rather, for 
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Turing machines, the involution would likely be a transformation on the tape (a CA 
involution?) along with a change in the current state of the machine. Finally, the 
locality of the time-reversing involution is not completely granted either: even in 
CA, it would be interesting to see what happens if that requirement is removed. 
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